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Abstract 

We theoretically study the dependence of the quasiparticle (QP) scattering rate r on the uniaxial anisotropy of a 
Fermi surface with changing the magnetic field angle qm- We consider the QP scattering due to the non-magnetic 
impurities inside a single vortex core. The field-angle dependence of the quasiparticle scattering rate -T(qm) is sensitive 
to the sign-change of the pair potential. We show that with increasing the two dimensionality of the system, /"(om) 
reflects more clearly whether there is the sign-change in the pair potential. 

Key words: Unconventional superconductors, Field-angle dependent measurement, Quasiparticle scattering rate, Phase 

sensitive probe 
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1. Introduction 

In anisotropic superconductors, the pair poten- 
tial has the anisotropy in the fc-space. For exam- 
ple, cuprate superconductors and the heavy fermion 
superconductor CeCoIn5 are experimentaly identi- 
fied as d 2; 2_j / 2-wave superconductors [1,2] . These su- 
perconductors have the zero of the pair potential 
(i.e. gap node) in 7r/4 direction measured from the 
a-axis. When the azimuthal angle goes across the 
gap node direction, the pair potential changes its 
sign. It is of great importance to elucidate the pair 
potential symmetry because it reflects the charac- 
teristics of the pairing interaction. Therefore, many 
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experiments have been carried out in order to in- 
vestigate the pairing symmetry. The field-angle re- 
solved thermal conductivity and specific heat mea- 
surments are powerful techniques, which can de- 
tect the anisotropy of the pair potential amplitude 
[2]. However, they are not phase-sensitive probe. In- 
formation on the sign change of the pair potential 
is crucial for clarifying the unconventional Cooper 
pairing (e.g., d^.^-wave pair). There are a few 
phase-sensitive probes [1,3], but an easier way to de- 
tect the sign change of the pair potential is required. 
In a previous paper, we reported that the field- 
angle («m) dependence of the flux-flow resistivity 
is phase-sensitive in the case of an isotropic spheri- 
cal Fermi surface (FS) [4]. In this paper, the depen- 
dence of the quasiparticle (QP) scattering rate r on 
the uniaxial anisotropy of a FS is theoretically in- 
vestigated under rotating the applied magnetic field 
H . We take into account the electronic structure re- 
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fleeting the crystal structure anisotropy by means of 
changing the dimensionality of the system. Then we 
consider a spheroidal FS, which is closer to realistic 
systems than a spherical FS (see Fig. ??). As a result 
of our calculation, we show that with increasing the 
uniaxial anisotropy of the system, r depends more 
strongly on om and P(«m) reflects more clearly the 
Cooper pairing symmetry including the sign-change 
of the pair potential. 

2. Formulation 




l k x k y 



We consider a single vortex core in which non- 
magnetic impurities are distributed randomly. r 
is obtained by calculating the imaginary part of 
the impurity-induced self energy S on the basis of 
the quasiclassical theory of superconductivity [6]. 
S appears as an energy shift of the vortex bound 
state in the denominator of the regular quasiclas- 
sical Green's function. The quasiclassical Green's 
functions are obtained analytically in the vicin- 
ity of a vortex core by using the Kramer-Pesch 
approximation [7]. We treat the impurity scatter- 
ing by using the t-matrix approximation [8]. We 
assume an s-wave scattering potential and con- 
sider the Born limit. The pair potential is given as 
A(r,fe F ) = f(r)A d(k F )e i4 ', where /(r) denotes 
the spatial variation of the pair potential, Ao is 
its bulk amplitude and (f> is the azimuthal angle. 
f(\r\ — > oo) = 1 and f(r = 0) = 0. The center of 
the vortex line is located at r = 0. d(k F ) represents 
the anisotropy of the pair potential in the fc-space. 
fep is the Fermi wave vector (i.e. the position on 
a spheroidal FS). The scattering rate r(e) for the 
vortex bound states with energy e is given by [9] 



£{s) 



(l-sgn[d(fc F )d(feF)] co s0) 

1 \VF±(K)\ \d(k F )\ 

\sme\\v F± (k F )\\d(k F )\ 



e -u(so,k F ) e -u(so,kp) 



FS'/ FS 



(1) 



where _T n is the scattering rate in the normal state. 
The brackets (• • • )fs an d (• ■ • )fs' mean the aver- 
ages on a FS with respect to k F and k F respec- 
tively, like (• ■ • ) FS = (l/iso)JdS F (k F )/\v F (k F )\ ■ ■ ■ 
with dS F being an area element on a FS and 
v o — J dS F (k F ) /\v F (k F )\. The QP is scattered from 



Fig. 1. Cross section view of a spheroidal FS along a certain 
rfik with the anisotropy ■y. The vertical axis is fez-axis and 
the horizontal axis is k x or fey-axis. The position on the 
spheroidal FS is characterized by the parameter (</)k ,0k) on 
a spherical shell with the radius fcp. 

the initial position fcp to the final one fc F on the 
spheroidal FS. |vFJ_(fcF)| and I^f^^f)! reflect the 
uniaxial anisotropy of the FS. The index g_Lhmeans 
the component projected onto the plane normal 
to H . The coherence length in the bulk is defined 
by Co = v F j_/irA with v F± = (\v F± (k F )\) FS - The 
details of Eq. (1) are discussed in Rcf. [4,9]. 

In this study, we consider the case where H is ro- 
tated in the a — b plane perpendicular to the k z -a,xis. 
The relations between the coordinates fixed to H 
(vortex coordinates) and those fixed to the crystal 
axes (crystal coordinates) are described in detail in 
the case of spherical FS in Ref. [4] . In a two dimen- 
sional system, when H is applied within the a — b 
plane, a distorted vortex core is formed. However, 
we neglect the vortex core distortion here and as- 
sume an isotropic vortex core. We set the variation 
of the pair potential in the plane normal to H, like 
f(r) = tanh(|r|/£ ). 

We consider a spheroidal FS uniaxially extended 
to the fc z -direction (see Fig. 2). For a spheroidal FS, 
the area element is given by 



dS F = fcpy cos 2 8 k + 7 2 sin 2 9 k sin. k d4> k d6 k ., (2) 

where the uniaxial anisotropy parameter is defined 
as 7 = \/m c /m a b- The mass m a b and m c charac- 
terize the charge transport within the a — b plane 
and in the c-axis, respectively. k F is the radius of a 
spherical shell (see Fig. 2). We assume the system is 
isotropic in the a — b plane. <f> k (9 k ) is the azimuthal 
(polar) angle on the spherical shell. The position on 
the spheroidal FS k F is uniquely identified by the 
position on this spherical shell (<j> k , 9 k ) (Fig. 2). 

On the spheroidal FS, the Fermi velocity is ex- 
pressed by [10] 



^fC^f) = VF(acos4>ks'm9k+bsm(j)ksm9k+c(l/^)cos9k), 

where a, b, c are orthogonal unit vectors character- 
izing the crtystal axes, v F is the Fermi velocity of in 
the a — b plane. After transforming the coordinate 
system from the crystal coordinates into the vortex 
coordinates, we can get the relations 
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v F 
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v F 
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(4) 
cos^fc, (5) 



sm9 v (kp) 



l' F 



aMJsrne/fe, 



(6) 



\v F± (k F )\ 

where k F = k F {o> cos 0^ sin 9k + b sin <f>k sin (?& + 
c."f cos 9 k ). 9 v (k F ) is the angle of the QP trajec- 
tory measured from the aM-axis. The aM-axis is 
unit vector spanning the plane normal to H . Here, 
H = \H\(a cos ocm + bsinaM)- All quantities are 
represented in the coordinate system fixed to the 
crystal axes. Therefore we can calculate the «m de- 
pendence of r in the same way as in the case of a 
spherical FS [4]. 

3. Results 

We consider the following three pair potential 
models on the spheroidal FS. (i) Line-node s-wave: 
d(k F ) = | cos 24>k sin 2 9k \ ■ (ii) Point-node s-wave 
[5]: d(k F ) = (1 + cos 4</>fc sin 0fc)/2. (iii) d-wave: 
d(k F ) = cos20fesin 9k- These pairing functions 
have the gap nodes on the spheroidal FS in the di- 
rections of 4> k = (1 + 2n)7r/4 with n = 0, 1, 2, 3. The 
anti-node directions correspond to (f>k — nn/A. In 
Figs. 2 - 4, we show the «m dependence of r. Each 
plot corresponds to the different 7. In each pairing 
state, the left-side panel shows the numerical result 
for the QP energy e = 0.1 Ao and the right-side one 
shows the result for e = 0.2A . 

As shown in Fig. 2, in the case of the line-node 
s-wave pair, r exhibits its minimum when H is ori- 
ented parallel to the gap-node direction (qm = f/4). 
With increasing 7, the minimum of r gets smaller. 
This behavior stands out in the higher energy. In an 
isotropic spherical FS, the oscillation amplitude of 
r is of the order 6% (7%) of the maximum value for 
£ = 0.1A (s = 0.2A ) [4]. On the other hand, in 
the spheroidal FS with 7 = 4, the amplitude is of 
the order 30% (50%). 

While the similar behavior is seen also in the case 
of the point-node s-wave pair, the difference appears 
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Fig. 2. Plot of the QP scattering rate r vs. the applied field 
angle om in the case of the line-node s-wave pair for (a) 
the QP energy e = O.lAo and (b) e = 0.2Ao- Each curve 
is plotted for different anisotropy 7. The vertical axis is 
normalized by maximum value for each curve. 

in the energy dependence of r(c«M)- In the line-node 
s-wave case, the minimum of r(ayi) gets smaller 
with increasing the QP energy e. However in the 
point-node s-wave case, it gets larger with increas- 
ing e. The other difference is seen in the oscillation 
pattern. In the case of the line-node s-wave pair, 
the cusp like structure appears as shown in Fig. 2, 
but in the case of the point-node s-wave, the cusp 
like structure disappears and the broad minimum of 
r(aij) appears (see Fig. 3). 

The similar calculations are conducted also for 
the d-wave pair (Fig. 4). The maximum of _T ap- 
pears when H is parallel to the node direction. With 
increasing 7, the peak of r becomes higher. In an 
isotropic system, the peak height is of the order 
120% (140%) of the minimum value [4]. On the other 
hand, in the spheroidal FS with 7 = 4, the peak 
height is of the order 150% (180%). 

In summary, the difference in _T(«m) between 
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Fig. 3. Plot of the QP scattering rate r vs. the applied field 
angle on in the case of the point-node s-wave pair for (a) 
the QP energy e = O.lAo and (b) e = 0.2Ao- Each curve 
is plotted for different anisotropy 7. The vertical axis is 
normalized by maximum value for each curve. 

the s-wave and the d-wave pair is also seen for the 
spheroidal FS. With increasing the two dimension- 
ality of the Fermi surface, the difference becomes 
prominent. 



4. Discussion and Conclusion 

We clarify that -T(o!m) reflects a sign-change of 
the pair potential more clearly in a spheroidal FS 
than in an isotropic spherical FS. With increasing 7, 
the difference of the behavior of r(au) between the 
s-wave and the d-wave pair becomes larger. There- 
fore, -T(a]yi) is phase-sensitive also for a spheroidal 
FS. Information on the uniaxial anisotropy of FS is 
included in vp±(kp) in Eq. (1). The different behav- 
ior of r(«M) between a spherical and a spheroidal 
FS comes from the different magnitudes of vp±(kp) 
and VF±(k' F ). 



Fig. 4. Plot of the QP scattering rate r vs. the applied 
field angle qm in the case of the d-wave pair for (a) the QP 
energy e = 0.1 Ao and (b) e = 0.2Ao- Each curve is plotted 
for different anisotropy 7. The vertical axis is normalized by 
minimum value for each curve. 

In moderately clean regime, the flux-flow resistiv- 
ity pf(T) is proportional to r(e = k-sT) [11]. One 
can obtain the information of P(«m) by measuring 
Pi with rotating H within the a — b plane. However, 
the strongly two dimensional system is not suitable 
for the measurment that we assume here, since if H 
is applied within the a — b plane, a Josephson vor- 
tex is formed [12], which is beyond our theoretical 
framework. In this work, we consider an isotropic 
vortex core. Elucidation of the effect of a vortex core 
distortion on r is left for a future problem. 
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